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1. Introduction 

Magnetic moment is a fundamental property of particles. It determines the dynamical response 
of a bound system to a soft external stimulus, and provides valuable insight into internal strong in- 
teraction structure. Efforts to compute the magnetic moment on the lattice come in two categories. 
One is the form factor method which involves three-point functions |[l], ^, ||, 0, ^ ^ The other 
is the background field method using only two-point functions (mass shifts) ||8|, ^ [l^ [Tl|]. The 
form factor method requires an extrapolation to zero momentum transfer Gm{Q^ = 0) due to the 
non- vanishing minimum discrete momentum on the lattice [|l^]. The background field method, on 
the other hand, accesses the magnetic moment directly and cleanly but is limited to static prop- 
erties due to the use of a static field. Here we report a calculation of the vector meson magnetic 



moments in this method, in parallel to a recent calculation in the form factor method [13]. It is 
an extension of our earlier work on baryon magnetic moments [ p^ ] and electric [15] and magnetic 
polarizabilities [ p^ ] in the same method. 

2. IMethod 

For a Dirac particle of spin s in uniform fields, 

E±=m±piB (2.1) 

where the upper sign means spin up and the lower sign means spin-down relative to the magnetic 
field, and ji = g^^s. We use the following method to extract the g factors, 

2m+m_ (£■+ - m+) - (£"_ -m_) 



m+ + m_ / eBs 

In order to place a magnetic field on the lattice, we construct an analogy to the continuum case. 
The fermion action is modified by the minimal coupling prescription 

D^=d^+gG^+qA^ (2.3) 

where q is the charge of the fermion field and is the vector potential describing the background 
field. On the lattice, the prescription amounts to multiplying a U(l) phase factor to the gauge links. 
Choosing Ay = Bx, a constant magnetic field B can be introduced in the z-direction. Then the phase 
factor is in the y-links 

Uy exp {iqa^Bx)Uy. (2.4) 

The computational demand can be divided into three categories. The first is & fully -dynamical 
calculation. For each value of external field, a new dynamical ensemble is needed that couples 
to u-quark (q=l/3), d-and s-quark (q=-2/3). This requires a Monte Carlo algorithm that can treat 
the three flavors separately. Quark propagators are then computed on the ensembles with matching 
values. This has not been attempted. The second can be termed re-weighting in which a perturbative 
expansion of action in terms of external field is performed (see Ref. [^] for a calculation of the 
neutron electric polarizability in this method). The third is U(l) quenched. No field is applied 
in the Monte-Carlo generation of the gauge fields, only in the valence quark propagation in the 
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given gauge background. In this case, any gauge ensemble can be used to compute valence quark 
propagators. 

We use standard Wilson actions on 24^ lattice at jS = 6.0, both SU(3) and U(l) quenched, 
and six kappa values K"=0.1515, 0.1525, 0.1535, 0.1540, 0.1545, 0.1555, corresponding to pion 
mass of 1015, 908, 794, 732, 667, 522 MeV. The critical value of kappa is K^=0.1571. The strange 
quark mass is set at K"=0.1535. The source location for the quark propagators is (x,y,z,t)=(12,l,l,2). 
We analyzed 87 configurations. The following five dimensionless numbers Tj = gBa^=+0.00036, 
-0.00072, +0.00144, -0.00288, -1-0.00576 give four small B fields (two positive, two negative) at 
eBa2=-0.00108, -1-0.00216, -0.00432, -1-0.00864 for both u and d (or s) quarks. These field values 
do not obey the quantization condition for periodicity since they cause too large a distortion to the 
system. To minimize the boundary effects, we work with fixed (or Dirichlet) b.c. in the x-direction 
and large Nx, so that quarks originating in the middle of the lattice have little chance of propagating 
to the edge. To eliminate the contamination from the even-power terms, we calculate mass shifts 
both in the field B and its reverse —B for each value of B, then take the difference and divide by 2. 
Another benefit of repeating the calculation with the field reversed is that by taking the average of 
5m{B) and 5m{—B) in the same dataset, one can eliminate the odd-powered terms in the mass shift. 



The coefficient of the leading quadratic term is directly related to the magnetic polarizability [16]. 



3. Interpolating field 

For the p+ meson, we use the polarized forms 

The interaction energies E± are extracted from the correlation functions 

{n±^±) = \ [{^xr]x) ± i {{fjxriy) - {fjyrijc)) + (%%)] • (3-2) 



Eq. ( |3.2| ) implies that the polarization comes from the imaginary parts of the off-diagonal correla- 
tion between x and y components. These imaginary parts are zero in the absence of the external 
field. Other vector mesons have similar forms with different quark content = iid, = uu — dd, 
^ = ss, K*^ = su, K* = lis, K*^ = ds — sd. By symmetry, the magnetic moments are expected to 
be related by 

Mp- = -j"p+, Mpo =0; (3.3) 

and 

fe*-=-fe*+, small. (3.4) 

These relations are borne out numerically in our calculations. 

4. Results and discussion 

Fig. [I] displays a typical effective mass plot for both the mass and the mass shifts. Good 
plateaus exist for all six quark masses. Our results are extracted from the time window 10 to 13. 
Fig. ^ shows the mass shifts, defined as 5 = g{eBs) from Eq. ( p^ , as a function of the field for the 
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Figure 1 

magnetic 



; Effective mass plot for the p+ mass (left) at zero field and mass shifts (right) at the weakest 
field in lattice units, corresponding to the heaviest and lightest quark masses. 
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Figure 2: Mass shifts for the p+ meson as a function of the magnetic field in lattice units at the six quark 
masses (lightest in the lower right corner). The slope of the mass shift at each quark mass gives the g factor 
corresponding to that quark mass. The line is a fit using only the two smallest B values. 



4 



Magnetic Moment of Vector Mesons 



Frank X. Lee 



p+ meson. The slope gives the g-factor. There is good linear behavior going through the origin for 
all the fields when the quark mass is heavy, an indication that contamination from the higher-power 
terms has been effectively eliminated by the {5{B) — 5{—B))/2 procedure. This is also confirmed 
numerically by the smallness of intercept as shown in the same figure. At the lightest quark mass 
(lower right comer), there is a slight deviation from linear behavior at the stronger fields. For this 
reason, we only use the two smallest field values to do the linear fit at all the quark masses. 

Fig. H shows the g-factors for the vector mesons as a function of pion mass squared. The lines 
are simple chiral fits using the ansatzs 

g = aQ+aimjc, (4.1) 

and 

g = ao+aim,^ + a2ml. (4.2) 

They serve to show that there is onset of non-analytic behavior as pion mass is lowered, so a linear 
extrapolation is probably not a good idea. But overall the g-factors have a fairly weak quark mass 
dependence. At large quark masses, the g-factor of p+ approaches 2, consistent with a previous 
lattice calculation using charge-overlap method Jls]]. Our results for p+ are slightly higher than 
those from the form factor method (see Fig. 8 in [|l^]). The results confirmed that gp- = —gp+ and 
g}^,- = —gK*+- We confirmed gpo = numerically (not shown). The results also show that as far 
as g-factors are concerned the p mesons are quite similar to their strange counterparts K* mesons. 

Note that the extracted g-factors are in the particle's natural magnetons. To convert them into 
magnetic moments in terms of the commonly-used nuclear magnetons (Hn), we need to scale the 
results by the factor 938/M where M is the mass of the particle measured in the same calculation at 
each quark mass. Fig. ^ shows the results for p+ and K*^. The different quark-mass dependence 
between p+ and K*^ mostly comes from that in their masses that are used to convert the g-factors to 
magnetic moments. The values at the chiral limit extrapolated from Eq. (22) are /ip+ = 3.25(3)/iAr 



and = 2.81(l)/XAr. There is no experimental information. Compared to the form factor method 
(see Fig.7 in [pj|]), our results are again a little higher. At the strange quark mass point (the 3rd data 
point from the left), the two coincide to give a prediction for the magnetic moment of the 0(1020) 
meson, /i0 = 2.07(7)/Ia?- 

Fig. ^ shows the results for K*^. Our results confirm the expectation that ;Ujf*o is small. It is 
positive when the d-quark is heavier than the s-quark, exactly zero when they are equal, and turns 
negative when the d-quark is lighter than the s-quark. The same behavior has been observed in the 
form factor method (see Fig. 1 1 in [IT: 



5. Conclusion 



In conclusion, we have computed the magnetic moment of vector mesons on the lattice using 
the background field method and standard lattice technology. Our results are consistent with those 
from the form factor method. There is no experimental information so the lattice results can serve 
as first-principles predictions. The calculation can be improved by providing a full account of 
systematic errors present in the results, such as finite-volume effects. In addition, there is a need to 
push the calculations to smaller pion masses so that reliable chiral extrapolations can be applied. 
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Figure 3: G-factors for the (left) and K* (right) mesons as a function of pion mass squared. The 2 lines 



are chiral fits according to Eq. (4.1 ) (dashed), Eq. (4.2) (dotted) 
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Figure 4: Magnetic moments (in nuclear magnetons) for p+ and K 



*+ 



Nonetheless, our results demonstrate that the method is robust and relatively cheap. Only mass 
shifts are required. This may facilitate the push to smaller pion masses, perhaps with the help 
of chiral fermions (overlap, domain-wall, twisted mass, ...). Finally, we await fully-dynamical 
background-field calculations in order to see the effects of the quenched approximation. 



Acknowledgment 

This work is supported in part by U.S. Department of Energy under grant DE-FG02-95ER40907. 
W.W. acknowledges a research leave from Baylor University. The computing resources at NERSC 
and JLab have been used. 



6 



Magnetic Moment of Vector Mesons 



Frank X. Lee 



o 



ao 
a 

a 
a 



o 
> 



0.15 



0.10 - 



o 0.05 



0.00 



-0.05 




Figure 5: Magnetic moments (in nuclear magnetons) for K 
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